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^_^ ABSTRACT 

'^ ' Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for 

r^ , models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain 

O |i in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an 

I ■ unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the 

O ' submillimeter < z < 0.2 local luminosity function and infrared bolometric (8-1000 fim) local luminosity density based on SPIRE 

-^ , data from the HerMES Herschel Key Program over 14.7 deg^. Flux measurements in the three SPIRE channels at 250, 350 and 

C/j , 500 /im are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy 

distribution of SPIRE sources and use the 1/Vmax estimator to provide the first constraints on the monochromatic 250, 350 and 

500 fj,m as well as on the infrared bolometric (8-1000 fira) local luminosity function based on Herschel data. We compare our results 

fvj , with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. 

►^ . Our measurement of the infrared bolometric (8-1000 ^m) local luminosity function suggests a flat slope at low luminosity, and the 

r~^ ' inferred local luminosity density, 1.311q'2i ■ 10* LqMpc"^, is consistent with the range of values reported in recent literature. 
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"iJ 1. Introduction both extended our reach to higher redshifts and improved 
^^ our constraints on the IR bolometric (8-1000 fim) luminos- 
2 Local lummosity functions (LLFs) are essential for the m- ^^^ ^^ detected sources. Due to the moderate volume per 
T: : terpretation of the evolution of cosmic sources at any wave- ^^-^ ^^^ ^^^^ ^-^^^ one can probe at low redshifts, an essen- 
length. While the determination of high-redshift LFs re- y^j pre-requisite for determining the LLF is the imaging of 
qmres deep observations, LLFs can exploit shallower wider- j^^g^ g^j^^ ^l^i^ ^^^ ^^ ^^^ ^^^^ ^^gj^^l^ ^^^ ground-based 
, area sky maps. In addition, the extremely careful multi- submillimeter (SMM) cameras, and for example Dunne et 
^ , wavelength identification work required for high-redshift ^j^ (2OOO) carried out SCUBA observations of IRAS bright 
■ ■ . sources is greatly simplified by the lower areal density of galaxies instead. It was only recently that BLAST was able 
bright local sources. The pioneering exploration of the IR ^^ ^^^p j^^g^ ^^^^^ ^^ g^^ wavelengths and thus to first de- 
sky by the IRAS satellite revealed the rapid evolution of the termine the evolution of the SMM LLF (Eales et al. 2009) . 
infrared (IR) LLF (Saunders et al. 1990), 1 lustratmg the ^he much improved mapping speed and angular resolution 
importance of local studies at infrared (IR) wavelengths, ^f ^^e SPIRE (Griffin et al. 2010) instrument onboard the 
ISO was later able to follow the evolution of the IR LF jj^^.^hd Space Observatory (Pilbratt et al. 2010) is ex- 
up to z ~ 1 (Pozzi et al. 2004), while Spttzer studies us- ^^^-^^^^ ^y the Herschel Muhi-tiered Extragalactic Survey 
mg the MIPS 24 fim (Le Floc'h et al. 2005 Marleau et (HerMES, http://hermes.sussex.ac.uk) Key Program 
al. 2007, Rodighiero et al. 2010) and 70 fim (Prayer et al. ^q^-^^^ ^^ ^1_ 2010, in prep), which wiU eventually cover 
2006, Huynh et al. 2007, Magnelh et al. 2009) channels have jq ^^^2 ^j^^jg ^^-^^^ g^^jy Herschel studies are probing 

the evolution of SMM number counts and LFs up to high 

* Herschel is an ESA space observatory with science instru- redshifts (Oliver et al. 2010, Eales et al. 2010), we present a 

ments provided by European-led Principal Investigator consor- ^^.^^ assessment of the LLF as measured in the three SPIRE 
tia and with important participation from NASA. 



M. Vaccari et al.: The HerMES SPIRE submillimeter local luminosity function 



bands at 250, 350 and 500 /im, providing the first robust 
constraints on the Local Luminosity Density (LLD) at these 
wavelengths. This study will thus establish a reliable local 
benchmark at a very early stage of the Herschel mission 
against which models and high-redshift studies can be com- 
pared. Throughout, we adopt a standard cosmology with 
flM = 0.28, f^A = 0.72, Ho = 72 kms-^Mpc-i. 



2. Sample selection 

We use a small portion of the observations making up 
the HerMES Key Program, namely those taken dur- 
ing the Herschel Science Demonstration Phase (SDP) in 
September-October 2009 (Oliver et al. 2010, in prep). We 
employ an Herschel source catalog obtained following a 
technique newly developed within the HerMES consortium. 
To minimize the effect of source blending, the SPIRE fluxes 
are estimated via linear inversion methods, using the posi- 
tions of known 24 /im sources as a prior. The full method 
and resulting catalogs are described in Roseboom et al. 
(2010, in prep). The area jointly covered by Spitzer and 
Herschel is of 10.7 deg^ within the Lockman Hole (LH) and 
4.0 deg^ within the Spitzer Extragalactic First Look Survey 
(XFLS), a total area of 14.7 deg^. We select S/N > 5 
reliable SPIRE detections as defined by Roseboom et al. 
(2010, in prep). While sources not appearing in the input 
24 /im source list will thus not appear in our sample, this 
incompleteness is negligible at < z < 0.2 and the adopted 
SPIRE flux limits for all but the most extreme SEDs, which 
we do not expect in significant numbers at low redshift. As 
discussed by Roseboom et al. 2010 (in prep), the 24 /im 
prior thus ensures that the SPIRE catalogs are effectively 
complete and do not suffer from severe blending issues over 
the < z < 0.2 redshift range down to the adopted flux cuts 
of 40 mJy and 30 mJy in all bands in LH and XFLS, respec- 
tively. For all Herschel sources, multi-wavelength archival 
photometry from the ultraviolet to the far-infrared and op- 
tical spectroscopy is provided by the Spitzer-selected multi- 
wavelength catalog by Vaccari et al. (2010, in prep) cover- 
ing the HerMES wide-area fields. Available spectroscopic 
redshifts are mostly from SDSS and NED, but are comple- 
mented by those obtained by Huang et al. (2010, in prep) 
and Sajina et al. (2010, in prep). We use spectroscopic red- 
shifts when available, and photometric redshifts from SDSS 
DR7 (Abazajian et al. 2009) computed following the neu- 
ral network technique by Oyaizu et al. (2008). Fig. [T] shows 
how optical counterparts to our SMM sources are almost 
invariably at r < 21, thus safely within the SDSS nominal 
magnitude limit of r ~ 22.4. 



3. SED fitting 

Because Herschel for the first time allows the sampling of 
the SMM part of the spectrum at multiple wavelengths, 
one can better constrain the nature of Herschel sources 
by fitting their optical-to-SMM spectral energy distribution 
(SED) with spectral templates. This allows us not only to 
determine the required iiT-correction but also to constrain 
the IR bolometric luminosity of individual sources. For this 
work, we use SDSS ugriz, 2MASS JHKs, IRAC 4-band 
and MIPS 3-band photometry which are available over the 
whole 14.7 deg^. We modeled (see Fig.[2|) the 250 /im sam- 
ple using the hyperz code of Bolzonella et al. (2000) and 
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Fig. 1. Optical magnitude vs 250 /im flux for the full sam- 
ple. The SDSS r ~ 22.4 magnitude limit is indicated by the 
red line. 
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Fig. 2. Example of optical-to-SMM SED fitting. The de- 
tailed SED fitting allows us to determine the X-correction 
which is then applied on the basis of the best-fit SEDs 
and the IR bolometric (8-1000 /im) luminosity for the full 
250 /im sample. We used SDSS, 2MASS, IRAC, MIPS and 
SPIRE measurements. 



the SED templates of Polletta et al. (2007), including the 
slightly modified versions of a few of these templates in- 
troduced by Gruppioni et al. (2010). The latter authors 
account for the first indications from their SED fitting of 
Herschel sources by introducing a higher FIR bump which 
improves the fit to their PACS data points. Keeping the 
redshift of any given source fixed, we thus determined the 
iiT-correction at 250, 350 and 500 /im on the basis of the 
best-fit SEDs and the IR bolometric (8-1000 /im) luminos- 
ity for the full 250 /tm sample. Interestingly, we find that 
late-type SED templates such as Sd and Sdm provide the 
best fit to the majority of sources in our sample. However, 
a detailed analysis of the LLF as a function of SED type is 
beyond the scope of this work. 



4. LLF estimation 

In order to determine the monochromatic LLF at 250, 
350 and 500 /im, we employed the standard 1/Vmax LF 
estimator (Schmidt 1968). Each source is assigned the 
maximum redshift z^ax at which it would still appear 
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"within our flux-limited sample and the associated Vmax = 
■min{V{z < Zmax),V{z < 0.2)). For every luminosity bin 
the LLF estimate is cj) — J^i ^^ — ^^^ its associated vari- 



^i Vm 



ance is a% 



E. 



< 1 — ^^j y^ — . Using a similar approach, we then 

computed the IR Bolometric (8-1000 /xm) LLF integrating 
the best-fit SED for each source. For each luminosity bin 
we then calculated our best LLF estimate and associated 
error as the weighted average and the associated weighted 
error of measurements in the two fields. Our estimates of 
the monochromatic and IR bolometric LLFs are reported in 
Table m We performed a least-squares fit to these measure- 
ments with a modified Schechter (or double exponential) 
function as defined by Saunders et al. (1990) 
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to determine the IR Bolometric LLD contributed by SPIRE 
sources through its integration over the full LLF. Best-fit 
parameters are log C [Mpc~^] = —2.16, \og{Li,/LQ) = 
9.94, a = 1.00 and a = 0.50 and the estimated IR bolo- 



metric LLD is 1.31 



5. Results 
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The LLF predictions from recent literature and previous 
LLF estimates are compared with our monochromatic and 
IR bolometric LLF measurements in Fig. |3l The 250 and 
350 /im LLF measurements in the two fields agree at an 
average level of 15%, consistent with the levels of cosmic 
variance predicted by theoretical models for fields of this 
size at < z < 0.2 (Moster et al. 2010). They also agree 
reasonably well with the BLAST measurements of Eales et 
al. (2009) but have much greater precision and fewer possi- 
ble systematic problems. Measurements at 500 /im, based 
on a much smaller sample, show a significantly larger vari- 
ance between the two fields and thus should be treated with 
caution. Our measurements are slightly above most model- 
ing predictions we considered at all wavelengths, and par- 
ticularly at the luminosities close to the knee of the LLF 
which we probed with greater statistical significance. At 
these luminosities, however, Xu et al. (2001) better repro- 
duce the observed space densities than Negrello et al. (2007) 
and Franceschini et al. (2010). Our measurements of the 
IR LLF suggest a flatter low-luminosity slope than previ- 
ously reported by either Sanders et al. (2003) for an IRAS 
60-^m-selected sample or Rodighiero et al. (2010) for a 24- 
/im-selected one, but agrees with predictions by Lagache et 
al. (2004). However, this flatter slope is mostly driven by 
our lowest luminosity bin, where the sample size is small 
and we may start to be affected by incompleteness. Our 
estimate of 1.31^0:2^ ■ 10* ^oMpc-^ for the IR LLD at an 
average redshift of z ~ 0.1 excellently agrees with the val- 
ues obtained by Sanders et al. (2003), Magnelli et al. (2009) 
and Rodighiero et al. (2010), but is higher than those ob- 
tained by Le Floc'h et al. (2005) and Seymour et al. (2010). 
Previous measurements at both the low and high end of this 
range in LLD used different combinations of MIPS 24 /im, 
IRAS 60 /j,m and MIPS 70 /im selection. While we can now 
get a better handle on the whole of the SED thanks to using 
Herschel as well as Spitzer measurements, the small sample 
size at both the faint and bright end does not allow us to 
place stronger constraints on the IR LLD as yet. 



Table 1. Number of < z < 0.2 sources used for the 
determination of LLFs. The number of sources with (spec- 
troscopic/photometric) redshifts is indicated after the total 
number of sources. The 250 /J,m sample is also used to de- 
termine the IR bolometric LLF. Field areas are in deg^. 



Field 250 ^m 



500 /iin Area 



350 /im 

116 (99/17) __ ,__,_, 

94 (80/14) 27 (22/5) 4.0 



18 (16/2) 107 



TR 478 (369/109) 

XFLS 275 (222/53) 94 (80/14) 27 (22/5) 4.0 

Total 753 (591/162) 210 (179/31) 45 (38/7) 14.7 

Table 2. Monochromatic and IR bolometric LLFs (0 < 
z < 0.2). L indicates i' L^ for the monochromatic LLF and 
Lju for the IR bolometric LLF. L at the bin centers are 
expressed in Lq while LLF estimates and their errors are 
in [Mpc~^dex^^]. 

logL log(0, o-)250 log (<^, o-)350 log (0, (7)500 log(0, o-)/fl 
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-2.38 , 


-4.50 
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-3.21 , 


-3.95 


11.5 


- 




- 




- 




-4.24 , 


-4.65 



6. Conclusions 

Using SPIRE imaging from the HerMES Key Program cov- 
ering 14.7 deg^ (about 20% of the total SPIRE area even- 
tually to be covered by HerMES) we presented the flrst 
Herschel measurement of the SMM LLF. We computed the 
IR monochromatic 250, 350 and 500 /im as well as the IR 
bolometric (8-1000 /xm) LFs of SPIRE < z < 0.2 sources 
and compared them with modeling predictions for IR LLFs. 
Relying on Spitzer 24 /im positions to identify SMM sources 
and on optical spectroscopic surveys to determine their dis- 
tances, we placed solid statistical constraints on the SMM 
LLF and LLD at a very early stage of the Herschel mission 
and compared our results with pve-Herschel predictions. 
We found a slightly more abundant local submillimeter pop- 
ulation than predicted by a number of models, a flat slope 
of the IR LLF at low luminosity and an IR bolometric LLD 
of 1.31^0:2^ • lO^L0Mpc-3, toward the higher end of val- 
ues from recent literature. Once HerMES observations are 
completed, we will be able to better probe the faint and 
bright end of the SMM LLF and study its evolution with 
redshift in much greater detail. 
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Fig. 3. SPIRE 250, 350 and 500 [m\ and IR Bolometric LLF in LH+XFLS. In the monochromatic LLF panels, predictions 
by Xu et al. (2001) and Negrello et al. (2007) are indicated by dashed and dot-dot-dot-dashed lines respectively while 
model populations by Franceschini et al. (2010) are: dashed blue line are spirals, dot-dashed cyan line are starbursts, 
dashed red line are high-luminosity starbursts, dot-dot-dot-dashed green line are Type-I AGNs, and the black solid line 
is the total prediction. Spirals are predicted to make up the bulk of local sources detected by SPIRE, as confirmed by 
SED fitting. In the IR bolometric LLF panel, the red long-dashed line is the best-fitting modified Schechter (or double 
exponential) function to the IR bolometric LLF, diamonds and squares are measurements based on the z ~ IRAS 
60-/Lim-selected RBGS by Sanders et al. (2003) and on a < z < 0.3 Spitzer 24-/.tm-selected sample by Rodighiero et al. 
(2010). Also in the IR Bolometric LLF panel, black dot-dashed, dot-dot-dot-dashed and dotted lines indicate predictions 
by Lagache et al. (2004), Negrello et al. (2007) and Vahante et al. (2009). 
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